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SUMMARY 

This  report  considers  the  theoretical  analysis  of  measurement  errors  introduced  in  an 
underwater  stereopsis  system  which  has  front  window  misalignment.  It  develops  new 
range  and  size  equations  and  models  the  effects  of  window  misalignment,  range  equation 
approximation  and  system  sampling  effects.  A  scheme  is  given  which  will  allow 
underwater  calibration  in  a  practical  system.  The  report  concludes  that  the  new  range  and 
size  equations  developed  provide  enhanced  accuracy  to  a  point  where  the  system  sampling 
effects  become  the  dominant  source  of  error. 
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1.  INTRODUCTION 

The  pnnciple  of  using  a  pair  of  similar  TV  cameras  to  derive  the  range  and 
size  of  viewed  objects  is  well  known  I'  -f  This  report  analyses  and  describes  a 
problem  caused  by  front  window  misalignment  in  a  real  underwater  stereo  ranging 
system  and  presents  a  way  to  compensate  for  the  range  and  size  measurement  errors 
which  occur.  Whilst  previous  work  has  looked  at  compensating  for  optical  distortions 
such  schemes  are  often  impractical  in  a  real-time  system  and  do  not  consider  the 
effects  of  a  water  medium  in  front  of  the  camera.  One  paper  which  has  looked  at 
underwater  photogrammetry  1^'  does  not  examine  the  effects  upon  a  stereoscopic 
system,  nor  does  it  consider  the  effects  of  window  misalignment.  This  paper 
represents  a  pragmatic  approach  to  a  specific  problem  with  a  solution  which  is 
achievable  in  real-time  and  with  a  calibration  scheme  which  may  be  performed 
underwater. 

Figure  1  shows  a  typical  system  in  which  the  range  of  an  object  seen  by  both 
cameras  may  be  determined  using  triangulation.  Once  the  range  is  known  the  object 
size  may  be  determined  by  using  the  range  as  a  scaling  factor. 

With  an  air  medium  between  the  objects  and  cameras  the  range  and  size 
calculations  are  fairly  straightforward  but  when  the  medium  is  water  other  effects 
become  apparent.  Previous  experience  with  underwater  stereopsis  has  revealed  a 
problem  with  camera  alignment.  It  was  found  that  after  aligning  the  two  cameras  in 
air  (for  parallel  optical  axes)  they  became  optically  misaligned  when  submerged  It 
was  suspected  that  this  effect  was  due  to  a  misalignment  (or  tilt')  of  the  front  window 
of  the  cameras'  water-tight  housings  relative  to  the  cameras'  optical  axes.  Corena  and 
Gillyon  I’l  modelled  the  arrangement  by  ray  tracing  and  were  able  to  show  that  a 
constant  angular  offset  in  the  principal  ray  of  the  order  of  those  observed  (around 
0.3'’)  will  occur  with  front  window  tilts  of  around  1  “.  These  misalignments  may  occur 
during  manufacture,  assembly  or  during  use  and  there  is  usually  no  convenient  means 
of  correction. 

This  repon  considers  the  effect  of  these  misalignments  on  the  system  ranging 
performance.  First  the  normal  range  equation  is  derived  and  then  extended  to  take 
account  of  misalignment  errors.  Next,  the  way  in  which  an  underwater  calibration 
may  be  used  to  reduce  the  effects  of  misalignment  is  given.  Following  on  from  this, 
the  results  of  system  simulations  for  various  formulations  of  the  ranging  calculation 
are  presented  to  illustrate  how  the  gross  effects  of  misalignment  may  be  taken  into 
account  and  appropriate  compensation  made.  These  simulations  include  the  effects  of 
camera  focal  plane  sampling  to  assess  the  contribution  which  this  makes  to  the  range 
and  size  calculation  accuracy. 

2.  DERIVATION  OF  RANGE  AND  SIZE  FORMULAE 
2.1.  Simple  Range  and  Size  Formulae 

This  section  gives  the  derivation  of  the  simple  range  and  size  formulae  in 
which  no  account  is  taken  of  window  misalignment.  The  calculations  are  based 
around  an  optical  system  of  the  type  given  in  Figure  2a.  Note  that  anticlockwise 
angles  are  positive  and  the  coordinate  system  origin  lies  on  the  left  camera  optical  axis 
at  the  lens  centre. 

From  this  diagram  we  have 

B  =  0,-0, 
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where  B  is  the  cameras'  baseline 
length 

R  is  the  object's  range. 


So. 


{ana,  -  tana, 


R-.^  • 

</  ■ 


(1) 

(2) 

where  /  is  the  cameras'  focal  length 
i,  and  i,  are  the  left  and  right 
camera  image  positions  for 
the  object. 


In  practice  i,  and  i,  are  measured  from  the  TV  picture,  often  by  overlaying  the 
left  and  right  images  and  allowing  an  operator  to  'merge'  the  two  images  of  the  object 
being  ranged.  The  amount  of  shifting  of  the  images  relates  to  (i,  -  ij. 


m 


♦ 


To  derive  the  size  formula  we  begin  with  (from  Figure  2b) 

S  =  O,  -  O, 

where  S  is  the  object  size. 

=  RtanO,  -  Rtana, 

S  =  R(tana,  -  tana,) 


•  •  • 


Now. 

tana,  =  ^  and  likewise  for  tana, .  So  we  can  say  :• 

S=j(ij-i,)  (3) 

where  i,  and  i,  are  the  image 
positions  for  the  object 
sizing  points 

These  are  the  commonly  used  formulations  of  the  range  and  size  calculations. 
We  wiU  now  consider  the  effects  of  window  misalignment 

2.2.  Range  and  Size  Formulae  with  Window  Misalignment 

As  discussed  above,  underwater  systems  are  known  to  have  errors  introduced 
by  misaligiunent  of  the  front  window  of  the  water-tight  housing.  Appendix  A  provides 
a  theoretical  development  which  demonstrates  how  a  fixed  window  tilt  will  produce  a 
fixed  angular  offset  in  any  ray.  This  angular  offset  will  be  approximately  equal  to  one 
third  the  window  tilt  and  in  the  opposite  direction.  The  appendix  also  shows  how  the 
translational  shift  in  a  ray  will  be  negligible.  Thus  we  make  the  assumption  of  a  fixed 
tilt  error  during  our  development  of  the  new  range  and  size  equations. 


*  All  derivations  used  in  this  report  assume  a  camera  focused  at  infinity  and  therefore  use 
focal  length  and  image  distance  interchangeably.  For  a  typical  lens  of  8  mm  focal  length  this 
results  in  errors  of  0.08%  for  a  camera  focussed  at  10  m  and  0.8%  for  a  camera  focussed  at  1 
m. 
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Figure  3  shows  a  stereopsis  system  with  window  misalignment.  Note  that 
angular  offsets  are  positive  anticlockwise. 

Whilst  equation  (1)  holds  true,  equation  (2)  does  not.  since  the  image 
positions  are  no  longer  solely  dependent  upon  the  object  position  and  the  lens  focal 
length.  Hence  we  need  to  derive  the  range  and  size  equations  for  this  more  realistic 
system  model.  Appendix  B  gives  the  derivations  of  these  equations  which  are  simply 
stated  here.  For  the  range 


J  A  -A  ‘‘iAAr  '' 

A,  +  ‘At  +  —— — 

ii  -  U  -*•  /( A,  -  A, ) (t,  - 1, ) 


where  A,  and  A,  are  the  left  and 
right  camera  tilts  in  the 
horizontal  plane 


We  will  show  below  how  (A,  •  A,)  may  be  measured  by  performing  a 
calibration.  In  a  practical,  in-service  system  it  would  be  extremely  difficult  to  measure 
4  and  A,  independently  and  so  we  would  like  to  be  able  to  ignore  terms  in  which  they 
occur  independently  in  order  that  we  might  use  the  approximation  :• 


*  J 

in  which  all  terms  are  either  known  or  are  measurable.  We  will  show 
later  by  simulation  the  effects  of  ignoring  these  terms. 

Again  referring  to  appendix  B  and  Figure  4.  the  size  formula  is 

„  _ h’i,  * 


/+A/i>iV+ 


Since  the  A,  term  is  not  readily  measured  in  practice,  we  would  like  to 
approximate  this  by  equation  3  :• 

c  S.r  I 


We  will  later  show  by  simulation  the  effects  of  this  approximation  and  of  a 
further  approximation  of  the  range  equation; 


A) 
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3.  SYSTEM  CALIBRATION 

The  (X  ■  X)  term  in  the  range  equation  (equation  5)  may  be  measured  by 
performing  an  underwater  calibration.  A  target  with  crosshair  markers  on  the  same 
hon/ontal  line  and  separated  by  the  camera  baseline  distance.  B.  is  used.  Figure  5 
shows  such  a  target.  The  left  camera  is  boresighted  on  the  left  crosshair  and  the 
distance  from  boresight  of  the  nght  crosshair,  viewed  with  the  right  camera,  is  noted. 
This  IS  shown  m  Figure  6.  The  angle  (A,  -  X)  can  be  directly  measured  where 

I S,  ■  S,)  =  arcran  J  (S) 

1,  IS  measured  from  the  TV  picture/signal  w-ith  an 
operator  merging  images  as  described  above. 

Note  that  this  calibration  is  different  from  simply  applying  a  shift  to  the 
relative  image  offset.  The  -  X)  term  in  equation  S  is  a  shift  in  relative  image 

A) 

displacement  but  the - ~J - term  represents  a  shift  which  varies  with  image 

position. 

4.  SYSTEM  MODELLING 

System  modelling  was  performed  to  examine  the  bkely  accuracy  of  the  system 
using  the  approximated  range  and  size  equations  derived  earlier.  The  minimum 
required  accuracy  was  ±10%  at  10  m.  though  any  improvement  was  considered  highly 
beneficial.  For  all  calculations  a  Im  wide  target  at  10  m  range  was  assumed  to  be 
placed  at  250  positions  on  either  side  (on  a  horizontal  plane)  of  the  left  camera  optical 
axis.  This  corresponds  to  object  positions  throughout  the  cameras'  horizontal  fields  of 
view.  For  each  object  position  the  image  position  for  each  camera  was  calculated, 
taking  into  account  the  effects  of  the  front-window  misalignment  on  the  cameras. 
These  image  positions  were  then  used  in  the  various  range  and  size  equations  being 
modelled. 

Results  are  shown  graphically  in  Figures  7  to  1 1  with  the  object  position 
plotted  on  the  abscissa  and  range  or  size  ploned  on  the  ordinate.  In  all  cases  two 
examples  of  window  misalignment  are  given.  The  leftmost  plot  shows  the  case  where 
the  left  camera  angular  shift  in  the  horizontal  plane  is  -0.006  radians  (-0.34  degrees), 
the  right  plot  shows  the  case  where  the  left  camera  angular  shift  is  +0.006  radians.  In 
all  cases  the  right  camera  angular  shift  is  -^O.OOO  radians.  Hence  the  right  plot  is  the 
case  where  the  rays  are  shifted  equally  for  both  cameras. 

4.1.  Simple  Range  Equation 

Figure  7  shows  the  calculated  range  using  the  simple  range  equation  given  as 
equation  2  and  the  corresponding  calculated  size  using  equation  3  (which  uses  the 
calculated  range)  for  the  two  cases  of  window  misalignment  described  above. 

Bearing  in  mind  the  true  range  of  -10  m  it  can  be  seen  that  when  the  two 
cameras  are  not  similarly  misaligned  the  simple  range  formula  of  equation  2  is 
hopelessly  inaccurate  and  varies  throughout  the  camera  horizontal  field  of  view.  The 
effects  of  this  range  error  upon  the  size  calculation  are  also  apparent.  The  simple 
range  formula  does  not  meet  the  required  accuracy. 

When  the  cameras  are  similarly  misaligned  (right  side  plots)  the  simple  range 
formula  is  accurate  to  approximately  0.5%.  The  size  equation  can  be  seen  to  be  highly 


-4- 

UNCLASSIHED 


•  •  • 


•  •  • 


J 


UNCLASSIFIED 


SRL-0I30-TR 


accurate,  with  calculation  precision  effects  clearly  visible  on  the  amplified  scale  of  the 
plot. 

4.2.  Full  Range  and  Size  Formulae 

To  veiify  the  full  range  and  size  formulae  given  as  equations  4  and  6.  the 
same  simulations  were  performed  using  these  formulae.  The  results  are  shown  in 
Figure  8.  For  the  two  misalignment  conditions  descnbed  above  this  figure  shows 
calculated  range,  percentage  range  error,  calculated  size  and  percentage  size  error. 

The  figure  shows  the  range  and  size  error  to  be  no  worse  than  0  0006%  and 
almost  uniform  throughout  the  field  of  view.  This  error  may  be  due  to  the  small  angle 
approximation  used  in  the  equations'  denvation. 

The  complete  range  and  size  formulae  are  more  than  adequate  to  meet  the 
system  accuracy  requirements  though  as  mentioned  previously  several  terms  are  not 
readily  obtainable  in  a  practical  system. 

4.3.  Range  and  Size  Formula  Approximations 

To  assess  which  terms  of  the  full  range  formula  it  is  possible  to  drop,  the 
above  simulation  was  run  for  two  cut-down'  versions  of  the  range  formula  of  equation 
4.  These  were  equation  7: 


and  equation  S  ;• 


JL 


<r  A, -A,) 


Equation  7  represents  the  simple  range  formula  with  a  shift  term.  Equation  S 
represents  the  simple  range  formula  with  a  shift  term  and  an  image  position 
dependent  shift  term. 

Figure  9  shows  the  calculated  range  and  range  error  for  these  formulae  and 
for  the  full  range  formula  (equation  4)  as  well  as  size  and  size  error  using  equation  3. 
For  size  calculations  the  range  value  used  is  the  calculated  range  not  the  true  range. 
On  some  graphs  only  two  plots  are  apparent.  In  these  the  plots  corresponding  to 
equations  4  and  5  are  coincident  though  with  an  expanded  scale  they  would  not  be. 
with  the  plot  for  equation  5  having  a  similar  shape  to  that  for  equation  7. 

The  plots  for  equation  7  ate  the  least  accurate  with  accuracy  decreasing  with 
distance  from  the  camera  optical  axis.  All  equations  meet  the  accuracy  requirement 
though  equation  5  offers  significant  improvement  without  requiring  additional 
information. 


As  with  previous  results  all  equations  ate  more  accurate  when  window 
misalignments  for  the  two  cameras  ate  the  same. 

Having  shown  that  it  is  possible  :o  derive  range  and  size  to  the  required 
accuracy  using  practically  useable  approximations  to  the  true  formulae  we  now 
consider  the  effect  of  data  sampling  on  range  and  size  accuracy.  Since  equation  5 
offers  an  improvement  over  equation  7  at  no  expense  (except  for  slightly  more 
computation)  this  formula  will  be  used  in  the  further  analysis. 
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4.4.  Sampling  Effects 

This  section  identifies  the  likely  errors  ansing  from  the  use  of  a  sampled 
system,  le  a  CCD  TV  camera.  Simulations  were  performed  as  above  but  the  calculated 
image  position  o  each  object  position  and  each  camera  was  subjected  to  binning. 
The  bin  siz^  .d  was  the  dimension  of  a  single  CCD  detector  element.  A  calibrauon 
was  pert-  -^icd  by  similarly  binmng  the  (A,  -  values  for  the  various  tilts.  Following 
this  me  range  using  equation  S  and  the  size  using  equation  3  were  calculated. 

A  typical  CCD  camera  will  have  up  to  768  elements  honzontally  and 
Simulations  were  initially  performed  for  1600  object  positions  throughout  the  field  of 
view  to  avoid  subsampling  artifacts.  For  the  purposes  of  illustration  the  graphs 
presented  here  were  produced  using  500  object  positions  as  for  the  earlier  simulations 
since  individual  lines  are  not  visible  with  1600  plot  points.  The  plots  shown  are 
representative  of  those  obtained  with  more  samples. 

The  top  half  of  Figure  10  shows  range  and  range  error  for  the  two  standard 
lilts  as  used  above.  The  effects  of  sampling  are  clearly  apparent.  Notice  that  the 
fluctuations  in  calculated  range  are  grouped  together  and  that  a  transition  between 
groups  occurs  where  the  binmng  causes  the  image  position  error  to  jump'  a  whole 
pixel. 

For  the  case  where  the  left  camera  tilt  is  different  from  the  nght  camera  a 
distinctive  bending  can  be  observed.  The  sampling  has  in  fact  accentuated  the 
bending  in  the  baste  range  formula  -  an  effect  which  was  invisible  in  the  non-sampled 
case  plotted  in  Figure  d. 

The  error  in  range  is  -7.5%  to  >2.5%  which  is  within  the  system  accuracy 
requirements. 

The  lower  plots  in  Figure  10  show  how  halving  the  bin  size  improves  the 
range  error.  This  may  be  achieved  in  practice  by  interpolating  between  adjacent  image 
pixels  using  image  processing  techniques,  assuming  that  the  image  of  the  ranging 
point  spans  two  adjacent  pixels.  The  range  error  for  the  smaller  bin  size  is 
approximately  t2.5%  -  well  within  the  system  requirements. 

It  IS  interesting  to  note  how  in  a  non-sampled  system  and  where  the  left  and 
right  cameras  had  the  same  tilt  the  range  equation  (equation  5)  produced  a  symmetric 
plot.  In  a  sampled  system  however  the  grouping'  of  the  fluctuations  in  calculated 
range  is  asymmetric  and  vanes  with  the  alt. 

4.5.  Effect  of  Range  Calculation  on  Size  Measurement 

Throughout  this  section  results  of  size  measurement  have  been  presented  for 
various  range  formulae  though  all  have  used  the  same  size  formula,  given  as  equation 
3.  Whilst  it  is  clear  that  the  range  calculation  has  an  impact  upon  the  size  calculation 
no  quantitative  assessment  of  this  effect  has  been  made.  Presented  here  are  the  results 
of  two  simulations.  One  uses  the  range  calculation  of  equation  5  to  derive  range 
values  which  are  then  used  in  the  equation  3  size  formula.  The  other  uses  the  true  ( 10 
m)  range  directly  in  the  equation  3  size  calculation.  In  both  cases  sampling  effects  are 
included  with  the  bin  size  equal  to  half  the  CCD  detector  size  (ie  we  assume  some 
interpolation). 

Figure  1 1  gives  the  results  of  the  two  simulations  with  the  normal  two  tilt 
cases  for  each.  The  upper  plots  show  the  simulation  where  calculated  range  is  used  in 
the  size  formula.  The  lower  plots  u.se  the  true  (uncalculated  and  unsampled)  range. 

The  effect  of  the  range  calculation  upon  the  size  calculation  is  quite  clear. 
Both  the  charactenstic  curve  (for  left  tilt  «  right  tilt)  of  the  range  formula  and  the 
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lampling  effects  impose  themselves  upon  the  calculated  size,  as  might  be  expected. 
Clearly  a  funher  improvement  in  the  range  calculation  wol  J  be  reflected  in  the  size 
accuracy. 


5.  CONCLUSIONS 

This  report  has  shown  how.  in  a  practical  underwater  stereo  ranging  and 
sizing  system,  the  standard  range  and  size  formulae  may  not  be  sufficiently  accurate  to 
allow  useful  measurements. 

it  has  demonstrated  how  the  effects  of  window  misalignment  may  be  made 
insignificant  through  the  use  of  an  electronic  underwater  calibration  scheme  and  a 
new  ranging  formulation. 

The  effects  of  image  sampling  upon  range  and  size  accuracy  have  been 
demonstrated  and  now  represent  the  major  limitations  upon  range  and  size  accuracy 
even  when  image  pixel  interpolation  is  allowed.  It  may  be  possible  to  funher 
interpolate  the  stereo  images  to  allow  image  matching  to  one  quaner  of  a  pixel  (one 
quaner  of  a  CCD  element)  accuracy  but  this  will  require  further  investigation. 

Simulations  were  performed  at  a  a  typical  operating  range  of  10  m.  Ranges 
much  in  excess  of  this  are  only  achievable  in  very  clear  water.  For  such  ranges  the 
sampling  effects  desenbed  here  will  become  more  dominant  if  the  camera  baseline  is 
not  increased. 

The  theoretical  results  presented  here  are  currently  being  tested 
expenmentally. 
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Figure  1  •  Ranging  With  Stereo  TV  Cameras 
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Figure  5  -  Calibration  Chan 
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Figure  7  -  Calculated  Range  and  Size  for  the  Simple  Formula 
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Figure  8  -  Range  and  Size  Calculations  for  Full  Formula 
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Figure  9  -  Comparison  of  Range  Equations 


Range  and  size  plots  using  Equations  4.  5  A.?.  NOTE :  Plots  for  eqations  4  AS  are  sometimes 
coincident. 
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Figure  10  •  Range  Calculation  with  two  Sample  Sizes 
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Figure  11  •  Effect  of  Calculated  Range  Inaccuracy  on  the  Size  Calculation  •  • 
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APPENDIX  A .  THEORETICAL  DEVELOPMENT  OF  THE  OPTICAL 
EFFECTS  OF  WINDOW  MISALIGNMENT  IN  A  THREE  MEDIA 

SYSTEM 


For  an  optical  system  in  water  we  will  show  how  the  deviation  in  any  light  ray 
caused  by  a  fixed  tilt  in  the  system  optics  front  window  produces  a  fixed  angular 
offset  in  perceived  angle  of  incidence.  We  also  show  that  the  shift  in  the  position  of 
the  emergent  ray  is  so  small  that  it  may  be  ignored. 

Figure  A1  shows  a  light  ray  travelling  through  three  media  which  in  our 
system  will  typically  be  water,  glass  (the  front  window)  and  air.  Angles  are  positive 
anticlockwise  and  a  conventional  cartesian  coordinate  system  is  assumed.  These 
media  have  refractive  indices  N|.  N2.  and  N3.  The  window  has  thickness  t  and  the 
angles  of  incidence  and  refraction  are  mariced  as  0.  This  system  has  no  window  tilt. 

System  without  window  tilt 

Angular  deviation  of  emergent  ray. 

N  N, 

Firstly.  let  k,  =  — ^  and  it,  =  — ^ 

N,  ■  N, 

Also,  let  T,  ~  tan6,  for  n  =  1,2,3. 

A  convenient  paraxial  approximation  to  Snell’s  law  gives;- 

and 

and  substituting  from  (A.l)  into  (A.2) 

T,  =  kM  - 


(A.1) 

(A.2) 

(A.3) 


•  • 


Deviation  in  emergent  ray  position,  Y3 

we  can  also  say  that 

Y,  =  Y^+tT,  -  (A.4) 

Substituting  from  (A.l)  into  (A.3) 

Y,  =  Y^+tkJ^  -  (A.5)  ^ 

Equations  (A.3)  and  (A.5)  give  the  angle  and  position  of  the  emergent  ray  for 
a  system  with  no  window  tilt.  We  will  now  consider  the  effects  of  window  tilt  on 
these  parameters. 


C 
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System  with  window  tilt 

Angular  deviation  of  emergent  ray. 

Figure  A2  shows  the  situation  with  a  tilted  window.  Primed  angles  (eg  6,  ) 
represent  measurements  relative  to  a  coordinate  system  which  has  rotated  with  the 

window.  Unprimed  angles  are  measured  relative  to  the  system  optical  axis.  K,  is  the 

new  value  of  K,.  due  to  the  window  tilt.  The  window  has  been  rotated  about  the 
position  where  the  incident  ray  impinges  on  the  front  of  the  window.  We  must 
establish  the  angular  deviation  in  the  emergent  ray  caused  by  the  window  tilt.  From 

the  figure.  6,'  is  : 

0,  =0, -A  -  (A.6) 

where  A  is  the  tilt  angle,  positive  anticlockwise. 

So 


7j  =  tan(9|  -  A) 

=  7j  -  A  in  the  paraxial  approximation  (A.7) 


Now.  from  (A.3) : 

Tj  =  k.kjT;  (A.8) 

From  Figure  A2  we  can  also  say  : 

r,  =s  T,  +fanA 
=  )k,kj7J  +ranA 
From  (A.7) : 

Ty  a  -  A)  -t-  A  in  the  paraxial  approximation 

=  k,kj7I -A(k,kj -1)  (A.9) 


For  a  water/glass/air  system  with  =  — 


Hence,  comparing  this  result  with  equation  (A.3).  the  deviation  in  any  ray 
caused  by  the  window  tilt  is  approximately  equal  to  1/3  the  angle  of  tilt  and  in  the 
opposite  direction.  This  concurs  with  the  experimentally  determined  value  obtained 
by  Corena  and  Gillyon.  Note  that  for  an  air/glass/air  system  there  is  no  angular 
deviation  caused  by  the  window  tilt.  This  also  corresponds  with  field  observations. 
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Deviation  in  emergent  ray  position,  Y3 
With  a  tilt,  referring  to  Figure  A2.  we  have : 

AB  =  ( 

AC  =  -^ 
cosQ, 

CD  =  ACsin(d2  +A) 
rsinjO,  +A) 
cosQ, 

_  f(sin6,  cos  A  +  cosQ,  sin  A  ) 
cosQ^ 

=  l(  tanQ^  cos  A  +  sin  A) 

t 

=  r(fan0j(l-5in' A>-  +sinA) 


In  the  paraxial  approximatioa  tan6  =>63:  sin8,  and  sin'  6^0,  so 


CD  =  f(r, +A) 

—  iT,  +’lA 

Hence 

y,  =  +tt  +fA 

Using  (A.l) : 

Kj  =  »;  +r/i:,7;  +tA 

Using  (A. 7): 

Yj=Y,+tk,{Tt-A)  +  tA 

Yj  =  ■¥ tk^T^+t Ail -k,}  (A.10) 


With  a  watei/glass/air  system  k^==  —  so,  comparing  (A.  10)  with  (A.3),  the 


At 

shift  in  Y  position  caused  by  the  tilted  window  will  be  approximately  —  ■  Using 

typical  values  of  ^  =  1*  and  r^O.OIm  the  shift  caused  by  the  tilted  window  is 
approximately  I9pm  and  therefore  we  ignore  it  in  our  modelling. 


‘3j 


if 


At: 


•  • 


•  • 


c 
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From  which  we  can  obtain: 

iiiX\  -  A) 

‘i  A)  + - J - ■  ',) 


The  formula  for  the  object  range. 

Derivation  of  Size  Formula 

From  figure  4,  the  object  size,  S  is 

5  =Oj  -  O, 

Using  the  same  form  of  derivation  as  for  equation  (3).  we  obtain: 
5  =  R(tana2  -  tana/} 


We  can  also  say  :- 
1 1 

j  =  tan((Xi  + 

and  using  the  derivations  for  tanOL,  and  rnna,  we  can  obtain  :- 


tana,  -  ■ 


and  similarly  for  tan  a. 


•  • 


Substituting  into  (B.3)  we  obtain 


'  L  LiAi  *  ii  *'/»' A  < 

f*  f  f  r  f  f 
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Simplifying  and  multiplying  cop  and  bottom  by  f 


f 


S  =  R 


t,-i,  -*•  Vfi,-i,l 
.  i,iA^ 

/+  ^ 


The  formula  for  the  object  size. 


(6) 
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